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We propose a novel quantum Monte Carlo method in configuration space, which stochastically 
samples the contribution from a large secondary space to the effective Hamiltonian in the energy 
dependent partitioning of Lowdin. The method treats quasi-degenerate electronic states on a target 
energy with bond dissociations and electronic excitations avoiding significant amount of the nega- 
tive sign problem. The performance is tested with small model systems of H4 and N2 at various 
configurations with quasi-degeneracy. 



I. INTRODUCTION 

Quantum Monte Carlo (QMC) methods in configura- 
tion space^ ^ have recently drawn attention of quantum 
chemists for obtaining accurate correlation energies of 
ah initio electronic structure theory. In these methods, 
the ground state wave functions are obtained stochasti- 
cally by a Monte Carlo simulation in the configuration 
space typically of Slater determinants. The auxiliary 
field QMC (AFQMC) method approximately controls the 
phase problem by writing the imaginary-time propagator 
of a many-body system with two-body interactions in 
terms of propagators for independent particles interact- 
ing with external AF.^^ More naive treatments are the 
projector Monte Carlo (PMC) method of Ohtsuka and 
Nagas^^M^ and the full configuration interaction QMC 
(FCIQMC) method of Booth, Alavi and coworkers,^ ^ 
which directly sample the coefficients in a FCI problem 
using walkers on configuration space basis functions. A 
long imaginary-time evolution of the wave function of 

^(r) = e-^(^-^)^(0) (1.1) 

shrinks the contributions of all states with finite excita- 
tion energies exponentially to magnify the ground state 
solution. Then the propagation for a short imaginary- 
time is calculated by a simulation of the differential equa- 
tion, 

^ = -(F-^)*(r). (1.2) 

Apart from the capability for applications, the schemes 
of the population dynamics of walkers in PMC and 
FCIQMC are similar except for the treatment of the di- 
agonal contribution of Hamiltonian matrix. 

The QMC methods suffer from a sign problem which 
prohibits stable energy convergence. The initiator vari- 
ant of FCIQMC^ (i-FCIQMC) overcomes the problem 
by controlling the growth of walkers of sign-incoherent 
progeny in an approximate manner with a bias of sur- 
viving criteria on a matrix with a sign problem. It is 
demonstrated that i-FCIQMC reduces the prefactor of 
exponential scaling significantly,*^ and the method has 



been successfully applied to the calculations of first-row 
diatomic molecules^ and even to those of real solids .^^ 
More recently, z-FCIQMC has been generalized in such a 
way that the propagation of a set of important determi- 
nants is performed by deterministic projection, resulting 
in a less severe sign problem and a large reduction in 
the statistical error. The basis set truncation error can 
be reduced by combining the QMC methods with F12 
electronic structure theory. A method for obtaining 
isolated excited states within the FCIQMC framework 
has been also proposed very recently.^ 

Despite the effort devoted to the advancements, none 
of the aforementioned QMC methods is effective for elec- 
tronic states in a quasi-degenerated or degenerated sit- 
uation. This is because those QMC deal with only a 
single-root problem in the Hilbert space. In this paper, 
we propose a novel QMC scheme to calculate the effec- 
tive Hamiltonian of the Lowdin partitioning rather than 
a single-root energy. The method treats multiple solu- 
tions near a target energy in a model space simultanously. 
This fact enables us to apply the new QMC method for 
various interesting phenomena in chemistry and physics 
involving bond dissociations and electronic excitations. 



II. METHOD 

A. Effective Hamiltonian and imaginary-time 
evolution of the transfer matrix 

In this section, we formulate the basic properties of 
effective Hamiltonian and its imaginary-time evolut ion. 
Let us begin with the Lowdin partitioning technique 
for the exact effective Hamiltonian. The entire configu- 
ration space is divided into the model space P contain- 
ing solution (s) of interest and its complement Q. The 
Schrodingier equation in the partitioned basis is in the 
form 
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For the state energy the CI vectors Cp and Cc 
exphcitly related as 



Cq = TqpCp, 
Tqp -(Hqq 
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where Tqp is the transfer matrix, and (2.5) is cahed the 
energy dependent partitioning (EDP) condition. Cq is 
then eUminated from the Schrodinger equation as 
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with the effective Hamiltonian, 
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According to the expansion of the model space, the cou- 
pling between P- and Q-spaces through the resolvent in 
(2.5) is diminished systematically as the gap between E 



and the eigen values of Hqq increases. The EDP tech- 
nique is applicable to any excited state by changing the 
target energy E as long as the principal component of 
the state is involved in the model space. In this case, Cp 
for different E are non-orthogonal except for the perfect 
degeneracy. 

One notes that an effective Hamiltonian appointing all 
solutions in the P-space universally is determined by the 
eigenvalue independ ent p artitioning (EIP) with the same 
number of variables 

Hqp + HqqTqp — Tqp(Hpp + HpqTqp) = 0. (2.9) 

We do not take this strategy in this paper because the 



numerical treatment of the linear equation (2.5) is more 
straightforward for the purpose of QMC implementation. 
The solution of the EIP equation is not unique due to 
the nonlinear structure of the form. If all of the states 
in the model space are well-separated from those in the 
Q-space, the transfer matrix Tqp can be expanded in a 
Ray leigh- Schrodinger perturbation series, which often 
accompanies intruder state problems. 

The imaginary-time evolution ( |1.2| ) consists of P- and 
Q-space contributions, ^ = ^p + ^q. The evolution 
of ^Q is obtained by multiplying the projector onto the 
Q-space, Q = 1 — P, 



dr 



-Q[{H - E)^ci{r) + H^p{t)], (2.10) 



or in the CI coefficients, 

-(Hqq - ^lQQ)CQ(r) - HQpCp(r). (2.11) 



^Cq 



dr 



Substituting (2.4) into the above, and breaking up the 



summation over Cp(r), we obtain the propagation of the 
transfer matrix. 



^TQp(r) 
dr 



-(Hqq - ^lQQ)TQp(r) - Hqp. (2.12) 



The stationary solution of th e eq uation dTQp{r)/dr = 
satisfies the EDP condition (|2.5|). Then the exact effec- 



tive Hamiltonian Hp?, is obtained from (2.7) and (2.8). 



For a given state energy, the TQp(r) propagation is de- 
coupled for each P-space determinant, and t he Np equa- 
tions can be treated independently. Eq. ( 2.12 ) bears 
close resemblance with the coefficients form of ( |1.2| ) ex- 
cept for the additional contribution of —Hqp. This fact 
enables us to develop a QMC algorithm in a manner sim- 
ilar to the other QMC methods in configuration space as 
explained below. 



B. Quantum Monte Carlo algorithm 

Henceforth, we use the notations /, J, ... and A, 5,... 
for the Slater determinants in the P- and Q-spaces, re- 
spectively. The dimension of the secondary space Q can 
be prohibitively large. We therefore intend to treat the 
imaginary-time evolution of TQp(r) stochastically with 
a small storage for Hp|> in the model space QMC algo- 
rithm (MSQMC). For this purpose, we augment (2.12) 
by multiplying a diagonal matrix Np^p containing fixed 
values of booster weights for the P-space determinants, 
and propagate it with a small time interval (5r, 



,5T'Qp(r) = 5r[p\^^{r) + p'^^{r) + p\^>{r)], (2.13) 

where Pq^\ Pqp diagonal (D), off-diagonal 

Q-space (Q), and P-space (P) contributions, respectively. 
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-H^"^T'Qp(r), 
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(2.14) 
(2.15) 
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and Hqq are diagonal and off-diagonal matrices 



of Hqq, and rQp(r) = TQp(r)N(];^. Although the 
booster weights can be fractional, we assume them pos- 
itive integers for simplicity in this work. The MSQMC 
method for the effective Hamiltonian represents T'Qp(r) 
as collections of walkers. 



(2.17) 



Sn and hhjj^ are the sign and corresponding determinant of 
the walker /i, respectively, and the distribution of walkers 
is dependent on r. The Monte Carlo steps correspond- 
ing to (2.14), (2.15), and (2.16) are independent for each 



walker as well as the P-space determinant /. Thus the 
merged walkers of different imaginary time steps repre- 
sent the solution of linear equation with much smaller 
number of walkers compared to the true dimension of 
TQp(r) as 



T'ai 



(2.18) 
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where ri and Tf are the initial and final times for the imag- 
inary time integration, and Ar = Tf — ri is the interval. 
For a fixed Nj^\ the number of walkers Nj^^ fluctuates 
during the imaginary- time evolution. The average num- 
ber of walkers TVj^^ is dependent on the coupling of / 
and the Q-space. Instantaneous Tai is obtained by re- 
moving the booster weie; ht, Tai = T'^i/Nf'\ We then 
accumulate Hpp with the instantaneous Tqp. In what 
follows is the practical algorithm of MSQMC, in which 
the individual steps of death/cloning, spawning, and an- 
nihilation are entirely parallel to those in FCIQMC.^ 

1. Secular part: Diagonalize Hp|> for Cp and E. 
Hp|> = Hpp for the flrst iteration. 

2. QMC part: All operations are reflected in a new 
list of walkers, and there is no interference until 
the annihilation step. 

(a) diagonal death/cloning step of pQ^^(r): Each 
walker in the Q-space dies with the probability 



St{H,. - E). 



(2.19) 



Although the cloning event for p^^^ < does 
not happen with a normal choice of a 
walker may be cloned with the probability 
The treatment of > 1 follows the 
way in the FCIQMC algorithm.^ 

(b) Q-space spawning step of p^\r): For each 
walker /i in the Q-space, a candidate is chosen 
randomly from the interacting determinants 
with the probability where is the 

number of Slater determinants in the Q-space 
interacting with hz excluding tz itself. Then a 
new walker u is spawned with the probability 



= -5rn.Hx 



(2.20) 



for pi^^ > 0. If pi^^ < 0, a walker with oppo- 
site sign is spawned with probability \p[^^ \ . 

(c) P-space spawning step of pQp(r): Inside a 
loop over a discretized index of the booster 
weight for each /, a determinant is chosen ran- 
domly as in the Q-space spawning, and a new 
walker in the Q-space u is spawned with the 
probability 



-SrniH^^i, 



(2.21) 



CP) 

for pu > 0, where n/ is the number of Slater 
determinants in the Q-space interacting with 
I. If pu < 0, a walker with opposite sign is 
spawned with probability \pu ^ | . 

(d) Annihilation step: Remove pair determinants 
with opposite signs from the list of walkers. 



(e) Sigma step: Using the list of new walkers, ac- 
cumulate Epp as 



1 



ArN 



(b) 



(2.22) 



for effective Ham iltonian Hp|) in the symmet- 
ric form of (2.7). 



For a probability p > 1, the event is caused int(p) 
times immediately with an additional event with 
the probability p—mt{p)W The steps (a) - (e) are 
repeated until the limit of the micro cycle. 

The secular and QMC parts are repeated during the 
macro iteration cycle. 



C. Stochastic promotion/ demotion stage 

For partitioning of the spaces in the MSQMC method, 
the P-space is preferred to be compact to shrink the com- 
putational cost. Contrarily, the presence of high-wighted 
amplitudes on determinants in the Q-space deteriorates 
the accuracy of MSQMC. In the static mode of MSQMC, 
the partitioning can be performed by selecting predomi- 
nating determinants from a complete active space (CAS) 
CI, which is generally unavailable for large molecules. 

For more general applications, we introduce a stochas- 
tic promotion/demotion (SPD) stage for the model space 
exploiting a sampling of QMC on top of a long production 
run. After certain QMC steps, SPD moves a determinant 
in the Q-space to the P-space if the number of walkers 
on the determinant exceeds TpA^'^*^ where Tp is a promo- 
tion threshold, and TV^*) = Ar(b) +7V(w) ig the sum of the 
booster weight and number of walkers. Nevertheless, the 
stochastic promotion can choose unwanted Slater deter- 
minants with small amplitudes. Thus, the determinants 
with amplitudes smaller than a threshold, < Td, are 
demoted in the subsequent secular step. 



III. RESULTS AND DISCUSSIONS 

We examine the performance of MSQMC on the model 
systems of H4 and N2 molecules, in which a full diagonal- 
ization is possible with a Hamiltonian matrix on mem- 
ory. SPD is performed with the contracted basis of the 
solution of Hpp rather than the individual Slater de- 
terminants using the flrst macro iteration. The lowest 
10 Slater determinants in Hjj are employed for the ini- 
tial Hpp of SPD. An N'^^^ constant simulation is per- 
formed during SPD. Population control is accomplished 
by the simple scaling = Ar(b)Ar(t)/(Ar(b) + A^W) 

every 10 steps. We use the default parameters of the 
promotion/demotion thresholds, Tp = 1 x 10~^ and 
Td = 1 X 10~^. The energy obtained from a diagonal- 
ization of Hpp is usually sufficiently low, and no pop- 
ulation control is employed for MSQMC after SPD in 
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FIG. 1: Behavior of the deviation of E from FCI for the H4 
model. Three simulations with different seeds are performed 
for each method and geometry. The SPD stage promoted 44 
and 50 Slater determinants for a = 0.0 and 0.5 respectively. 



this work. The total number of walkers is temporarily 
increased to A^^*^ = 2 x 10^ for SPD, while the booster 
weight is kept N^^^ = 5 x 10^ for all / in MSQMC. The 
death/cloning step is associated with the instantaneous 
energy until the integration is turned on at the 5th macro 
iteration. The definition of the instantaneous energy in 
this work is the expectation value of the instantaneous 
contribution of the distribution of walkers Eq. (2.17) to 
the effective Hamiltonian using CI coefficients fixed in 
each macro iteration. Each macro iteration contains 10^ 
steps with Sr = 0.01. Three simulations using different 
seeds for random numbers are performed for each method 
to confirm the magnitude of statistical error. We do not 
employ any bias in the sampling space of walkers to avoid 
introducing additional error albeit the initiator extension 
of MSQMC is straightforward 



TABLE I: Statistical measures {mEh) of all energies at every 
macro iteration cycle after 5 x 10^ QMC steps. 



FCIQMC 



MSQMC 



ATp = 1 ATp = 2 SPD 



a 

Aabs 
Amax 

Astd 

A (s) a 
^std 



0.5 



0.0 



a 

Aabs 

Amax 

Astd 



0.09 


0.19 




0.02 


0.19 


0.49 




0.06 


0.06 


0.12 




0.02 


0.12 


0.22 




0.03 


0.53 


0.36 


0.14 


0.01 


1.31 


0.77 


0.26 


0.03 


0.62 


0.37 


0.15 


0.02 


0.97 


0.38 


0.22 


0.01 



^Standard deviation of energies with different seeds at the 50th 
macro iteration step. 
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FIG. 2: Growth of the numbers of walkers in the first 1,000 
steps of the MSQMC/SPD simulation with 44 promoted 
Slater determinants in the model space for a = 0.0. 



A. H4 model 

The first example is the H4 model with the double zeta 
plus polarization basis set. This model has been em- 
ployed freque ntly to test multi-reference electronic struc- 
ture methods. 1^2"^ The non-dynamic correlation effects 
are significant near the regular square form with the angle 
parameter a = 0. The dimensions of the FCI problems 
are 5,068 and 5,524 for a = and 0.5, respectively, in 
D2h point group symmetry, and 10,036 for other a in 
C2V symmetry. Fig. [l] shows the behaviors of errors in 
E compared to FCI with respect to the number of QMC 



steps for a = 0.5 and 0.0. The number of walkers grows 
to nearly 4,000 for a = in the MSQMC simulation only 
with the Hartree-Fock (HF) determinant in the P-space 
(MSQMC/1). Correspondingly, we set the initial num- 
ber of walkers N^"^^ = 4, 000 in a FCIQMC simulation 
with the logarithmic control of energy shift. ^ A perfect 

iV(w) 

constant simulation is not possible with the poor 
initial HF energv, and A/'(^) finally grows by ca. 20% . 
For FCIQMC, the integrated value of local energy with 
intermediate normalization is used for comparison of the 
results. For a = 0.5, the HF determinant is predominat- 
ing, and all results are accurate to 0.5 mEi^ after 3 x 10^ 
steps. FCIQMC compares with MSQMC only with the 



5 




-♦-FCIQMC 


-□-MSQMC/1 


-ii-MSQMC/SPD 


^MRCI 


=-<-CCSD 


^3-MRexpT 


SS-MRCC 







0.2 

a 



FIG. 3: Errors of QMC and correlated methods for different 
a of the H4 model. The results for the deterministic methods 
are taken from Refs.^^^. 



HF determinant in the P-space denoted by MSQMC/1 
while MSQMC/SPD outperforms them significantly. Ta- 
ble [T| shows statistical measures of energy, mean absolute 
error Aabs, maximum absolute error Amax, and standard 
deviation Agtd. Agtd is in agreement with Aabs in each 

(s) 

simulation. A^^^ of energies with different seeds at the 
50th macro iteration step also compares with Agtd as 
a good candidate for the measure of error estimation. 
a = 0.0 is the most difficult geometry for H4. Clearly 
the accuracies of FCIQMC and MSQMC/1 are deterio- 
rated compared to the results for a = 0.5. In this ge- 



ometry, (a^)^(62n)^ and (a^)^(&3n)^ determinants in the 
D2h have equal weights. The extended model space with 
the 2 Slater determinants, MSQMC/2, clearly improves 
the result. MSQMC/SPD is as accurate as the result for 
a = 0.5. Based on Aabs and Agtd, MSQMC/SPD is 1- 
2 orders of magnitude more accurate than FCIQMC or 
MSQMC/1 especially for the quasi-degenerate situation. 

Fig. [2] shows the growing behavior of the numbers of 
walkers in MSQMC/SPD for a = 0. Each line corre- 
sponds to the number for one of the Slater determinants 
in the model space promoted at the SPD stage. The 
number distribution grows very rapidly in the initial 200 
steps. There are two groups in the equilibrated distri- 
bution around A/"^^^ = 1,000 and 1,500, indicating the 
presence of the main component of the wave function and 
determinants interacting relatively strongly with those in 
the Q-space, respectively. In either case, the numbers of 
walkers needed for the description of the Q-space contri- 
bution are noticeably reduced on average from 4,000 in 
the MSQMC/1 case. The accuracy of MSQMC/QMC is 
supported by the reduced coupling between the P- and 
Q-spaces. 

It is tempting to compare the accuracy of the QMC 
methods with those of correlated methods available in 
ab initio electronic structure theory. Fig. [3] compare re- 
sults for various a parameters. The averaged energy after 
5 X 10^ steps of simulations with different seeds is used 
as the resulting value of a scheme for all QMC methods. 
Compared to CCSD, all multi-reference methods supplies 
well-balanced dependency of the geometrical parameter 
a. The errors of FCIQMC and MSQMC/1 are smaller 
than those of all deterministic methods on average. Nev- 
ertheless, these QMC methods will not be extremely use- 
ful for practical applications due to the unsystematic er- 
rors. MSQMC/SPD significantly reduces the error; the 
result is accurate to a few 10 /i^h for all a. 



B. Bond stretching with electronic excitation of N2 



The second system is the N2 molecule at different bond 
distances. Despite the small dimension of the FCI prob- 
lem (8,152 in D2h symmetry), it is not easy for QMC 
to to deal with the system since several CI coefficients 
with equally large amplitudes are present at a stretched 
distance. We use the same parameters as the H4 case ex- 
cept for the MSQMC/1 simulation, in which the weight 
of HF determinant is insignificant at large r to increase 
N'^^h We therefore use the reduced value of N^^^ = 100 
at r = 4.2 and 6.0 ao to make the number of walkers 
comparable with those in FCIQMC. Such a significant 
growth does not occur in MSQMC/SPD. The FCI space 



is generated by CASSCF calculations with 10 electro ns in 
10 orbitals (10,10) with the aug-cc-pVDZ basis set.l^^ 
The orbitals are in the supplementary material.l^ 

The results for the ground and excited sate energies 
are summarized in Table [llj Although the exact solu- 
tions of the two lowlying states in the Ag symmetry are 
uncoupled and ^5]+, they cannot be separated by 
symmetry in the basis of Slater determinants. The num- 
ber of the model space determinants generated in the 
SPD stage ranges from 50 to 80 dependent on the dis- 
tance. Near the equilibrium geometry, r = 2.068 ao, the 
ground state is dominated by a single Slater determinant, 
and FCIQMC and MSQMC/1 reveal accurate. We have 
performed energy specific MSQMC simulations appoint- 
ing the ground and excited states separately as denoted 
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TABLE II: Ground and excited state energies of the N2 molecule at several bond distances. 



^Ag '2iA.g '2iAg \Ag 



r 2.068 ao 

FCI -109.142 767 -108.414 054 19.829 

FCIQMC 0.000 135 (0.000 185) 

MSQMC/1 0.000 287 (0.000 318) 

MSQMC/SPD/^i 0.000 022 (0.000 013) 0.016 824 (0.000 020) 0.457 

MSQMC/SPD/^2 -0.007 198 (0.000 161) -0.000 030 (0.000 096) 0.195 

MSQMC/SPD/^2 -^1" -0.001 

r - 4.2 ao 

FCI -108.807 006 -108.791 832 0.413 

FCIQMC 0.003 768 (0.003 758) 

MSQMC/1^ -0.001 082 (0.005 063) 

MSQMC/SPD/^i 0.000 083 (0.000 059) 0.002 117 (0.001 760) 0.055 

MSQMC/SPD/^2 -0.000 699 (0.000 430) 0.000 237 (0.000 281) 0.025 

MSQMC/SPD/^2 -^1" 0.004 

r — 6.0 ao 

FCI -108.801 215 -108.800 615 0.016 

FCIQMC 0.001 026 (0.004 316) 

MSQMC/1^ -0.006 170 (0.005 178) 

MSQMC/SPD/^i -0.000 076 (0.000 142) -0.000 073 (0.000 163) 0.000 

MSQMC/SPD/^2 -0.000 121 (0.000 210) 0.000 114 (0.000 127) 0.006 

MSQMC/SPD/^2 - ^1" 0.005 



^Excitation energy as the difference of the target energies of 
MSQMC/SPD/£;i and E2 simulations. 

^The reduced number of virtual walkers A^(^) = 100 is used to 
keep A^*^*^ comparable. 



by MSQMC/SPD/£;i and E2. The energy of the target 
state in each simulation is calculated precisely in spite of 
the fact that the SPD stage is performed for the ground 
state. The effective Hamiltonian is not universally accu- 
rate with somewhat large excitation energies due to the 
large energy gap between the states at r = 2.068 ao- The 
difference between the two energy specific results supplies 
the excitation energy very close to the FCI value. 



At the extended bond distance 4.2 ao, the energy gap 
is reduced to 0.4 eV. The quasi-degeneracy clearly makes 
FCIQMC and MSQMC/1 imprecise. A^'^j of the methods 
amount to several m£^h indicating a significant negative- 
sign problem induced from the neighboring states. The 
results of MSQMC/SPD retains similar accuracies for the 
target states. The reduction of the excitation energy in- 
creases the transferability of the effective Hamiltonian 
with the entire errors less than 0.1 eV. These features 
hold for the results at 6.0 ao- In this case, the energy gap 
is comparable to the static error of the MSQMC/SPD 
simulations. 



C. Perspectives on real applications 

The examination of the preliminary implementation 
of MSQMC motivated by the newly develope d form ula 
for the time evolution of the transfer matrix (2.12) has 



clearly demonstrated the effectiveness of the method 
for electronic states of small model systems with quasi- 
degeneracy and excitations. So the remaining question 
is, whether the effectiveness is preserved on scaling up 
to realistic systems with large CI dimensions of interest 
in chemistry and physics? Indeed, MSQMC, which only 
cures the instability from quasi-degeneracy by including 
interacting states into H^®^\ is not a sort of means to 
eliminate an exponential growth in the number of neces- 
sary walkers due to sign-incoherent spawning events from 
determinants with low populations for larger systems or 
larger basis sets. The initiator method introduced by the 
Alavi group^ ^ has ameliorated the situation in FCIQMC 
reducing the sign problem, which smoothly disappears as 
increasing the number of walkers. Likewise, it is expected 
that the initiator approach will complement MSQMC 
achieving the effectiveness for real systems. More specifi- 
cally, the initiator space at first equivalent to the P-space 
of MSQMC is enlarged such that a noninitiator determi- 
nant in the Q-space with a walker number exceeding a 
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threshold is made an initiator ehgible to contribute to 
the spawning processes onto unoccupied determinants. 
In this way, the initiator method wih be incorporated 
into MSQMC to produce an extremely useful scheme ca- 
pable of treating large and difficult systems with quasi- 
degeneracy and strong electron correlation. The initiator 
extension is meaningful for larger systems. A compari- 
son of the initiator variants, i-FCIQMC and i-MSQMC 
is therefore of interest and should be made with full-scale 
implementations . 

IV. CONCLUSIONS 

In this paper, we have proposed and tested the 
MSQMC method, which incorporates with the effective 
Hamiltonian of the Lowdin partitioning technique into a 
stochastic simulation for the FCI contribution outside the 
model space. The new method complements most of the 
defects in the previous QMC methods in configuration 
space. The main features of MSQMC are summarized in 
the following. 

1. The most severe negative-sign problem is originat- 
ing from a stochastic treatment of the competition 
of a small number of predominating Slater determi- 
nants. MSQMC sidesteps the problem by handling 
the principal component in a deterministic secular 
equation with an effective Hamiltonian in the same 
spirit as the semi-stochastic QMC method.^^^ 

2. MSQMC is capable of treating multi-root solutions 
in the vicinity of the target energy. This feature 
guarantees accurate descriptions of degenerate and 
quasi-degenerate electronic states involving a bond- 
breaking. 

3. The EDP condition holds for excited states with 
different target energies. MSQMC is capable of 
providing an accurate energy of excited state with a 
model space containing basis functions predominat- 
ing the state of interest. The effective Hamiltonian 
is transferable for states in a certain range of the 
target energy. 



4. The MSQMC steps are independent for each 
parental determinant in the model space. This 
fact enables an efficient parallel implementation 
by communicating the state energy and effective 
Hamiltonian in each macro iteration, in addition to 
the intrinsic scalability with respect to the number 
of processes with different seeds for random num- 
bers as well as the numbers of walkers in each of 
them. 

5. The stochastic promotion/demotion step is partic- 
ularly effective to extract important basis functions 
to describe the target state. 

MSQMC is computationally demanding compared to the 
standard electronic structure methods and other QMC 
methods for a single-root solution. Nevertheless, the re- 
quirement for a FCI problem shifts from the storage and 
communication of CI coefficients to computing resources, 
which can be easily obtained in modern computational 
environments. An implementation with on-t he-fly com- 
putation of Hamiltonian matrix elements for large-scale 
FCI problems will be reported elsewhere.— 
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